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^ . Abstract 

CL|| In Z decay into prompt charmonium, i.e. charmonium not originating 

Q^' from S-meson decays, the most important contribution is expected to 

O . come from colour-octet mechanisms. However, previous fixed-order 

calculations of the colour-octet contribution contain large logarithms 
which, in a more complete treatment, should be resummed to all or- 
^ ■ ders. We study this resummation by using a Monte Carlo QCD cas- 

^ I cade model and find that the fixed-order colour-octet result is dimin- 

ished by 15%. We compare the Monte Carlo calculations with results 
obtained by using analytical evolution equations. 



1 Introduction 



The production of charmonium and bottomonium states in various processes, 
especially at high-energy colliders [0, has recently received considerable ex- 
perimental and theoretical interest. New data have become available from 
pp p[, ep and e^e~ [||-|^ colliders. Theoretically, it has been realized 
that quarkonium production at colliders is dominated by parton fragmenta- 
tion [|^. For the factorization of the hard dynamics of QQ production and 
the soft dynamics of quarkonium bound state formation, a new framework 
is provided by non-relativistic QCD (NRQCD) P]. This new theoretical ap- 
proach has had initial successes and is now approaching a stage where its 
predictive power can really be judged by testing relations between various 
observables. 

Within NRQCD, quarkonium is produced either through colour- singlet 
mechanisms, where a perturbatively produced colour- singlet QQ pair evolves 
non-perturbatively into a quarkonium state, or through colour-octet mech- 
anisms, where the intermediate QQ pair is in a colour-octet state and the 
non-perturbative transition involves the exchange of soft gluons with the 
environment. Before the advent of NRQCD, quarkonium production was 
usually calculated taking into account only the colour-singlet channels. Con- 
tributions from colour-octet channels are usually suppressed by powers of t;^, 
the intrinsic velocity squared of the bound state, which apparently justifies 
the colour-singlet model. However, there are reactions where this suppression 
is either absent or compensated by an enhancement of colour-octet channels 
by powers of the strong coupling constant a^. 

A well-known example is the production of prompt ip' or direct prompt 
J /ip in pp collisions at large p^. (By "prompt" we mean those charmonium 
states which are not produced in 5-meson decays, and by "direct" those 
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which are not produced in the decays of other charmonium states.) The pro- 
duction rates measured by the CDF Collaboration are more than an order 
of magnitude above the colour- singlet-mo del prediction. This anomaly can 
be explained in the framework of NRQCD [T^, which includes a colour-octet 
production channel whose contribution is proportional to the probability of 
a non-perturbative transition, to be treated as a free parameter of the theory 
and fitted to the CDF data. 

The fact that NRQCD comes with a number of free parameters (although 
there are scaling rules which determine the order of magnitude of the transi- 
tion probabilities in terms of v"^) necessitates cross checks between different 



experiments. Colour-octet contributions to e.g. photoproduction []11| and 
fixed-target hadroproduction [|12[ have indeed been calculated, partly in or- 
der to provide such cross checks, partly as an attempt to explain discrepancies 



between the colour singlet model and experimental data [T^\. 

An interesting laboratory of quarkonium production is provided by e'^e~ 
annihilation at the pole. Experimental branching ratios or upper limits 
for decay into prompt quarkonium have recently been published by the 
OPAL 0,11 and DELPHI [|] Collaborations. Theoretical calculations 0,|14|- 
p!9| show that the largest contribution is expected to come from the same 
non-perturbative transition that in the NRQCD framework explains the ip' 
anomaly. Using the value of the NRQCD matrix element from Tevatron 
fits, the magnitude of this contribution is consistent with the experimental 



numbers |118|,|T9[ . However, the existing calculations have been done at a fixed 
order in ag, and the result includes large logarithms, ln(M|/Mjy^), which 
call for a resummation. It is the purpose of this paper to do this resummation 
and see how much the fixed-order results are changed. 

Our principal tool in studying the evolution effects is a Monte Carlo QCD 
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cascade model, implemented in the Ariadne program ^|. This approach 
was used in earlier work on the large p± hadroproduction of quarkonium. 

The evolution effects can also be studied using analytical evolution equa- 
tions p2[ to determine the scale dependence of fragmentation functions. 



Since the analytic fragmentation formalism fails in the threshold region where 
the colour-octet fragmentation process peaks, we derive our final results using 
the Monte Carlo method. Above the threshold region we use the analytical 
evolution equations as a check of the Monte Carlo implementation. 

In Section ^ below, we discuss in detail the colour- singlet and colour- 
octet mechanisms of parton fragmentation into charmonium and review the 
earlier work of other authors. We also describe the analytical fragmentation 
formalism used as a check of our Monte Carlo simulation. In Section ^ 
the Monte Carlo simulation of evolution effects with the Ariadne program 
is discussed. The Monte Carlo method is then applied to the colour-octet 
fragmentation process to obtain our main results. These are presented in 
Section ^, where we show that the evolution effects diminish the theoretical 
prediction of branching ratio into prompt charmonium by about 15%. 
Finally, the results are discussed in Section |^. 

2 QCD production mechanisms 

In decay into ^Si charmonium (which we denote generically as ip) the 
dominant colour- singlet contribution comes from charm fragmentation, i.e. 
the process where Z^ decay into cc is followed by the perturbative QCD pro- 
cess c — cc ^Si'^ + c and then the non-perturbative process cc ^S[^^ ijj. 
There is an equal contribution from c fragmentation. The cross section is 

^^Z^^ij + X) = 2 T{Z" ^ cc) Dfiz, /i), (1) 
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where z = E^/{Mz/2) is the fraction of the initial charm quark energy 
taken by the ip, Df{z,fi) is the charm fragmentation function into ip, and 



/i = 0{Mz/2) is the fragmentation scale [|l^]. The normalization of this 
contribution is related to the leptonic width of the ip. 



In [jT6[ large logarithms ln(M|/Mjy^) were resummed using Altarelli- 
Parisi evolution equations ||2^ for the scale dependence of Df{z,fi). Here 
we will use both Monte Carlo simulations and improved analytical evolution 
equations to study the evolution effects. 

The dominant colour-octet contribution comes from fragmentation mech- 
anisms where decay into qqg is followed by the perturbative splitting 



g cc 



3 0(8) 



Si and then the non-perturbative process cc 
leading contribution to the non-perturbative process comes from the emis- 
sion of two soft gluons. The probability of this non-perturbative transition 
is proportional to the NRQCD matrix element (0|(9g (^5*1) |0), i.e. the same 
element whose value has been fit to reproduce the CDF data. The normaliza- 
tion of the colour-octet component can therefore be obtained from Tevatron 
fits. 

In the following, all numbers were derived using the colour-octet matrix 
element {Q\Oi'^ {^Si) |0) = 0.0066 (GeV)^ obtained from a fit of the CDF 
data and the colour-sing let matrix element (0|Cf/'^(35i)|0) = 0.73 (GeV)^ 
TBI, related to the wavefunction at the origin of coordinate space -Rs(O) 



through |i?5(0)|2 = (27r/9)(0|O//'^(3^i)|0) = 0.512 (GeV)=^. Furthermore, 
rric = 1.5 GeV, = 2mc, and as{2mc) = 0.253 [|l^ were used. 

Below, we shall speak of decay into direct prompt J/ip- The produc- 
tion of Xc and ip' from cc ^S[^^ intermediate states and their decays into 
J/ip can be included by multiplying the octet contribution by a factor 2.1, 
derived from the octet matrix elements of Ref. and the branching ratios 
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given in |25|. The production of ifj' from cc '^S[^^ and its decay into J/ip can 



be included by multiplying the singlet contribution by a factor 1.3, derived 
using the colour-singlet-model relation 



\RisioW 



Mj/^j r( j/^ a) 



(2) 



with experimental leptonic widths 



This far, the calculation of the colour-octet contribution has only been 
done at a fixed order in perturbative QCD. The process is then represented 
by the Feynman diagram shown in Fig. |l]a and another diagram where the 
gluon is emitted by the antiquark line. The resulting charmonium energy 
spectrum is [JT^^ 
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which leads to an integrated width 
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giving Bt^^\Z^ direct prompt J/ip) = 0.68 ■ 10"^ for the J/ip. Hence this 
contribution is clearly larger than the colour- singlet charm fragmentation 
contribution Br^^\Z^ direct prompt J/ip) = 0.22 ■ 10"'^ |TB[. The z distri- 
bution is also markedly different, peaking around 2; = 0.1 whereas the singlet 
component peaks around z = 0.8. 
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Figure 1: The Feynman diagrams which describe the dominant colour-octet 
and colour- singlet mechanisms of fragmentation into prompt J/ip- (o-) 
Colour-octet gluon fragmentation. There is another diagram where the gluon 
is emitted by the antiquark line, (b) Colour-singlet charm fragmentation. 
There is (in an axial gauge) another diagram which represents c fragmenta- 
tion. 



Far from the threshold {z ^ M^/Mz) the differential decay width can be 
written as 



dV 

dz 
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The large logarithm arises from the nearly collinear splitting of an initial 
quark into a quark and a virtual gluon, the virtuality M| being small on the 
scale set by the Z^ mass. The double logarithm in Eq. @) arises from the 
soft gluon emission represented by the factors I/2; in Eq. (|^). In a complete 
calculation such logarithms should be resummed to all orders in perturbation 
theory. 

Just as in the colour singlet charm fragmentation case, such a resumma- 
tion can be performed by rewriting the decay width in terms of fragmentation 
functions and then resumming the large logarithms using analytical evolution 
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equations for the fragmentation functions. Such a treatment does, however, 
neglect terms of (9(M|/(z^M|)), which become important in the threshold 
region, where the colour octet J / ^/^-production peaks. We will therefore per- 
form the resummation by means of a Monte Carlo simulation, and use the 
analytical resummation described below only as check at large z. Note that 
even at the threshold the virtuality of the fragmenting quark can be as low 
as of the order M^Mz which is still much smaller than M|, thus motivating 
a fragmentation approach. 

In fragmentation language, the Feynman diagram in Fig. |T]a may be 
viewed either as quark production followed by quark fragmentation into J/il) 
(via splitting), or if the quark-gluon invariant mass is large, as gluon produc- 
tion followed by gluon fragmentation into J/ip. Correspondingly, the colour 
octet contribution can be written as p|, p!8| , |26[] 

dV , 
dz 



+ ^ {zx,fi)D]^{x,fi), 6 

JO X dZg ^ 

where /i sets the scale for the fragmentation functions and defines the minimal 
qg or qg invariant mass for the differential three-parton decay width 



dTjZ^ -^qqg) ,^ ^ ^ qq) 

dZa TC 







Inserting LO expressions for the fragmentation functions 
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in Eq. the LO large z result Eq. (j^) is reproduced up to corrections 
0(/.VM|). 

At scales /i of the order of M^, Eq. (||) is totally dominated by the quark 
fragmentation term. The corrections are 0{\\ but the leading logarithmic 
terms are still correctly reproduced. 

To resum the large logarithms we have used the improved evolution equa- 
tions [E21 



Dfiz,M'/z) = 0. (11) 



In the colour singlet case we have used the following charm fragmentation 
source term extracted from 



dt{z,xml) 



8a2|i?(0)|2 X 



(x^ - 2a; - 47) - z{x - l){x - 9) + 



4:z{l-z) , 4(8-7^-5^2) 



12^2 1-^ 



\2 



(12) 



(2 - 

Using this source term we have solved the evolution equation (|TD|) numerically 
and used Eq. (|I]) to find the colour singlet contribution to J/ip production 
in Z° decays. 

In the colour octet case the gluon fragmentation source term is a double 
delta function 

Equivalently, one may omit this source term and use Eq. (|^) as an initial 
condition at the scale 2mr. 
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We have solved the evolution equation (|TD|) numerically and used Eq. 
performing the convolution of the evolved gluon fragmentation function and 
the Z° — > qqg differential decay width Eq. (|^), to find the colour octet con- 
tribution to J/ip production in decays. 

The commonly used AP evolution differs from the improved evolution of 
Eq. (|1^) in two ways: 

• In the AP evolution equations the source term df^z,^"^) is put to zero 
and instead unphysical initial conditions are chosen at the scale 
such that the AP-evolved fragmentation functions agree with pertur- 
bative fragmentation functions at large scales. 

• In the AP evolution equations, the fragmentation functions on the right 
hand side of Eq. (|1^) are evaluated at the scale /x^ instead of at the 
scale yfi'^, and thus, fail to take into account the momentum constraint 
saying that the maximum virtuality of a parton taking a positive light- 
cone momentum fraction y from a parent parton with virtuality fi"^, is 

Neither of these points is crucial for our comparison with the Monte Carlo 
simulation at large scales and far from the threshold. 

3 The Monte Carlo method 

The implementation of charmonium production in the Ariadne event gen- 



erator was described in Ref. |21[], and here we will only give the main points 
and present some improvements. 

Three mechanisms are included (see below). In all cases, the perturbative 
splitting into a cc pair and the non-perturbative formation of the charmonium 
state are treated together as one step in the otherwise purely perturbative 
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QCD shower. The motivation behind including the soft formation of a char- 
monium state in the perturbative level, before the normal hadronization, is 
that all hard interactions in principle are already included in the splitting 
function. 

The three implemented mechanisms are treated somewhat differently: 

• Colour octet gluon fragmentation: implemented as a perturbative 
splitting of a gluon into a collinear colour octet cc ^S[^^ pair, followed 
by a non-perturbative transition into a ip and one soft gluon. To cor- 
rectly conserve quantum numbers, there should be at least two gluons 
emitted, but since the subsequent string fragmentation is stable w.r.t. 
addition of soft gluons, emitting only one to conserve colour should give 
a good description of the final state hadrons. For practical reasons, the 
delta function in the splitting kernel 6{z — 1) is given a small width by 
replacing it with step functions (1/e) 6{z — (1 — e)) 9{1 — z). It is now 
possible to study the final-state properties of the octet mechanism as 
a function of this width which typically is expected to be on the order 
of the squared velocity f ^ of the quarks inside the charmonium. 

• Colour singlet gluon fragmentation: a gluon is perturbatively split 
into a collinear colour sing let cc [^sf ^1 pair and two hard gluons. The 
cc pair directly forms a ip while the two gluons may continue cascad- 
ing until hadronization sets in. In Ref. [^] only one gluon was emit- 
ted to conserve colour and momentum, but now the full three-body 
phase space is explored with the two-gluon invariant mass correctly 
distributed according to Ref. |^. The final-state distributions dis- 
cussed in Ref. turned out to be insensitive to this change, and the 
conclusions therein still hold. The colour singlet gluon fragmentation 
contribution to J/ip production at LEP is very small and will not be 
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discussed further in this report. 

• Colour singlet charm quark fragmentation: a charm quark is split 
into a collinear cc '^Si^^ pair and a charm quark according to Fig. |l]b. 
The cc pair directly forms a ip while the charm quark may continue 
radiating. 

Besides the treatment of the soft gluons in the octet model, the main un- 
certainty in the implementation comes from the fact that the dipole cascade 
model in the Ariadne program in each step has all partons on-shell. In 
order for a gluon to split into a charmonium, some energy has to be borrowed 



from one of the colour-connected partons. As discussed in Ref. [21], this is a 
safe procedure as long as the cascade is strongly ordered, i.e. when the trans- 
verse momentum p±g of the gluon to be split is much larger than the trans- 
verse mass m±^ of the charmonium. For high-p^ charmonium production at 
hadron colliders this is a good approximation. But in decays there may 
even be situations where the mass m^+g of the formed charmonium-gluon 
state is larger than p_i_g, especially for the octet channel at small z. We have 
tried to estimate these uncertainties by e.g. explicitly requiring m^+g < p^g, 
and the results turned out to be insensitive to such modifications. There- 
fore we assume that the description is good also at small z, although further 
studies are needed. 

In this paper we will only use this program to study J/ip production at 
LEP, but the implementation is general enough for studying the fragmenta- 
tion production of any quarkonia at any experiment. 
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Figure 2: The z distribution of J /ip from colour- singlet charm quark frag- 
mentation. Solid and dashed histograms are from Ariadne in leading order 
and with full evolution respectively. The dashed line is from an analytical 
calculation to leading order, and the dotted and dash- dotted line are from 
analytical calculations with improved evolution using two different scales. 

4 Results 

In this section we present the results of our Monte Carlo simulations. We first 
compare with the analytical evolution, then we discuss the energy isolation 
cuts of the experimental analysis in the presence of evolution effects and 
finally proceed to determine the branching ratio of into prompt J/ip. 

In each case, we also discuss the dependence of our results on the en- 
ergy carried by the soft gluons emitted in the non-perturbative colour-octet 
transitions. This is of the order of the quark kinetic energies mcf^/2 in the 
charmonium rest frame and hence negligible in the first approximation. 

In Figs. and ^ we compare the results using the analytical fragmen- 
tation formalism with the results from the Monte Carlo simulation. The 
results agree well in the large z limit. Closer to the threshold region, where 
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Figure 3: The z distribution of J/ip from colour-octet gluon fragmentation. 
The solid line is from an analytical calculation to leading order. The dotted 
and dash-dotted lines are from analytical calculations with evolution using two 
different scales. The histograms are from Ariadne using different widths of 
the step function in the splitting function: long-dashed corresponds almost to 
a delta function (e = 10~^j, and short-dashed to a width of e = 0.2 ^ v"^ . 



the analytical fragmentation formalism is known to fail, the disagreement be- 
comes large, especially for octet gluon fragmentation. If we introduce a finite 
width in Eq. (0), corresponding to an energy of the order of M^v"^ carried 
by the soft gluons in the non-perturbative transition cc ^Sf'^ — ^ ip -\- gg, the 
distribution somewhat softened, as expected. 

In the experimental analysis, the separation of prompt J/ip from the B- 
decay component is primarily based on energy isolation cuts, supplemented 
by vertex detector information ||^, ^. The J/ip produced in 5-meson de- 
cays are typically accompanied by approximately coUinear, energetic light 
hadrons. The mechanisms of parton fragmentation into prompt J/ip, on the 
other hand, have a tendency to kick the cc pair rather far apart from the 
original parton direction, resulting in the production of charmonia with little 
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accompanying energy. Prompt charmonia can therefore be identified by re- 
quiring tlie accompanying energy -Econe within a cone of half-opening angle a 
around the charmonium direction to be less than some Eq. The LEP groups 
use ^0 = 4 GeV and a = 20° § or 30° [|. 

We have studied how much the evolution effects change the distribution 
of energy around the charmonium. The results are presented in Figs. ^ and 
^ where we plot the distribution of events as a function of accompanying 
energy, ((ir/(ii?cone)/r- In the majority of events, charmonium states com- 
ing from B decays are accompanied with energy -Econc > 4 GeV; the prompt 
component, on the other hand, contains a dominant peak at -Econe — 0, repre- 
senting events where there are no hard partons within 30° of the charmonium 
momentum direction. Although the fixed-order result for prompt production 
is changed by including the effects of the parton cascade, hadronization and 
a finite energy of the soft gluons emitted in the octet channel, the fixed-order 
calculation already gives a good idea of how the isolation cuts work in a re- 
alistic situation. Hence we do not expect any large corrections to the results 
of Refs. |] and [|. 

As an aside we want to point out that the threshold at z = 2M^/Mz is re- 
produced in the analytic fragmentation formalism if the distinction between 
the light cone momentum fraction y = p'^/Mz = {E^ + \p^\)/Mz and the 
energy fraction z = 2E^/Mz = y + y'^M^/M^ is taken into account. This 
distinction has been neglected in equations such as Eq. (|), since fragmenta- 
tion functions D{y,fi'^) are correct only up to corrections of O {y'"^ / M^) ■ 

Our results for the direct prompt J/ip energy spectrum dT/dz before and 
after the cut on Econe, are presented in Figs. |^ and 0. The shapes of the 
distributions are not significantly changed by the evolution effects. 

For the colour- singlet contribution from charm fragmentation (where evo- 
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Figure 4: The distribution of J/ ip events from charm fragmentation and from 
B decays as a function of the accompanying energy -Econe within a cone of 
half-opening angle a = 30° around the J/ip direction. Solid histogram is 
charm fragmentation from Ariadne to leading order on the parton level. 
Long- dashed histogram is the same but with full evolution and hadronization. 
Short-dashed histogram is J/ip from B decays. The distributions have been 
normalized to unity. 

lution effects can easily be included also in the analytical calculation by means 
of differential equations, see Fig. ||), the total integrated branching ratio is 
Br(i)(Z° direct prompt J ftp) = 0.27 ■ 10-^ which is is the same as the 
leading order result within the statistical error of 1% in the Monte Carlo 
calculation. This is consistent with the normal AP evolution which does not 
change the overall normalization. Also the improved evolution only differs 
from standard AP in the small-^ region, hence there is no large change in 
the normalization in that case either. Note, however, that for the charm 
fragmentation we have everywhere used a fixed as = 0.253 in the splitting 
function of Eq. (|T2|). By default, Ariadne would use a running also here, 
and with the Aqcd = 0.22 GeV fitted from global event shapes, the result is 
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Figure 5: The -Econe distribution as in Fig. ^. Solid histogram is gluon frag- 
mentation according to the octet model from Ariadne to leading order on 
the parton level. Dashed histograms are the same but with full evolution and 
hadronization: long-dashed with a delta function in the splitting function, and 
short-dashed a finite width of e = v"^. The dotted line is from an analytical 
calculation to leading order. 

20% higher. 

For the the colour-octet gluon fragmentation channel using Eq. ( [I3| ) we 
get Bt^^\Z^ —>■ direct prompt J/ip) = 0.58 ■ lO"'', which is 15% below the 
fixed-order analytical result of 0.68 ■ 10~^. Introducing a finite energy of 
the soft gluons in the octet mechanism lowers the value even further to 
Br^^\Z^ direct prompt J/ip) = 0.46 ■ 10~^. Note, however, that the 
value used for {0\O^^^ i^Si) |0) is from a fit ||2^ to Tevatron data neglecting 
the energy of the soft gluons. But high-p^ hadroproduction cross sections 
are approximately proportional to the fifth moment / dz z'^ D{z,fi) of the 
fragmentation function, so that a trigger bias is introduced and the normal- 
ization becomes very sensitive to the exact form of the fragmentation function 
at large z. It turns out that changing from a delta function in Eq. ([I3|) to 
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Figure 6: The z distribution of J/ip from charm fragmentation without (upper 
curves) and with (lower curves) isolation cut. Solid histograms are from 
Ariadne in leading order, and dashed histograms are from Ariadne with 
full evolution and hadronization. 

a finite width e = 0.2 ^ implies that (OlOg^^^ {^Si) |0) has to be doubled 
to give the same cross section at the Tevatron. This would mean that the 
branching ratio becomes Bt{Z^ — > direct prompt J/ip + X) = 0.92 ■ 10~^ in 
the finite width case. Also for the octet channel, the default running a.^ in 
the splitting function used in Ariadne gives a 20% higher branching ratio 
than the fixed as used here. 

To compare with experimental branching ratio of into total prompt 
J/ip, the colour-octet contribution has to be multiplied by a factor of 2.1 and 
the colour-singlet contribution by a factor of 1.3 to take into account Xc and 
tp' decays into J/ip. The resulting branching ratios are 1.57 ■ 10~^ for e = 
and 2.28 ■ 10~^ for e = 0.2, both consistent with the OPAL measurement ^ 
Br(Z° total prompt J/^ + X) = (1.9 ± 0.7 ± 0.5 ± 0.5) ■ 10"^ 
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Figure 7: The z distribution of J/ip from colour-octet gluon fragmentation 
without (upper curves) and with (lower curves) isolation cut. Solid lines are 
from an analytical leading order calculation. Dashed histograms are from 
Ariadne with full evolution and hadronization: long-dashed with a delta 
function in the splitting function, and short-dashed a finite width of e = v"^. 



5 Discussion 



We have studied evolution effects in Z° fragmentation into prompt J ftp by 
means of a Monte Carlo simulation. We worked within the framework of non- 
relativistic QCD (NRQCD), where the production of charmonium proceeds 
through either colour-octet or colour- singlet cc intermediate states. 

Evolution effects turn out to be relatively small, diminishing the colour- 
octet contribution to the decay width by 15% while colour- singlet contribu- 
tion is unchanged. The shape of the J/ip energy spectrum also turns out 
to be stable against the evolution effects. Thus the conclusion drawn by 



previous authors |T8|,|T9[ on the basis of a leading-order calculation - that 
the colour-octet contribution has larger normalization and a much softer en- 
ergy spectrum than the colour-singlet contribution - remains valid after the 
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inclusion of the evolution effects. 

In fact, the dominant uncertainties in the comparison between NRQCD 
predictions and experimentally measured branching ratios are due to uncer- 
tainties in the fits of NRQCD matrix elements and to the small statistics of 
the experimental results. 

The Tevatron fit of the {Q\0'l^^ i^Si) |0) matrix element is subject to a 
large uncertainty due to the trigger bias effect. The existing fits of Teva- 
tron data neglect the energy of the two soft gluons emitted in the transition 
cc ^Sf'' + gg, thus maximizing the hardness of the fragmentation 

function. A realistic emitted energy of the order of M^v"^ corresponds to a 
softening of the fragmentation function which suppresses the cross section 
by approximately a factor of two. The fit value of the matrix element would 
then have to be increased by this factor of two. 

Apart from this and other theoretical uncertainties, the statistical er- 
rors of the Tevatron data imply uncertainties in the matrix elements around 
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In the bottomonium sector, the approximation of a heavy quark mass 
is better than in the charmonium sector. On the other hand, there is less 
data available, and theoretical expressions involve more free parameters than 
in the charmonium sector. In the Tevatron experiments, contributions from 
Xb decays are not resolved from direct T production; in decays, even 
the different radial excitations of the T can currently not be resolved from 
each other. In principle, however, the analysis can be easily extended to the 



bottomonium sector 19 



The predictions of NRQCD are consistent with the observed quarkonium 
production in decays. However, both theoretical and experimental un- 
certainties are large. A precision test of NRQCD may become possible if 
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uncertainties in the Tevatron fits of NRQCD matrix elements are reduced 
and if more data become available at the pole. Evolution effects in the 
theoretical calculations are smaller than these uncertainties and are under 
control, as we have shown above. 
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